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The complete electroweak radiative 0{a) corrections to the Drell-Yan process at large invariant 
dimuon mass have been studied. All formulas for the cross sections and kinematical restrictions are 
presented in explicit form, for the simplification of calculation and coding the 9- and 5-functions 
are actively used. The FORTRAN code READY for the numerical analysis in the high energy 
region corresponding to the future experiments at the CERN Large Hadron Collider has been 
constructed. To simulate the detector acceptance we used the standard CMS detector cuts. The 
radiative corrections are found to become large at high dimuon mass M, the complete corrections at 
' "bare" setup change the dimuon mass distribution up to ~ —5.6% (—23.2%; —35.3%) at the LHC 

^ i energy and M = 1 (3; 5)TeV. 
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, ■ I. INTRODUCTION 



More than twenty years the Standard Model (SM) keeps for oneself the status of consistent and experimentally 
confirmed theory since the experimental data of past and present accelerators (LEP, SLC and Tevatron) have shown 
^sj ' no significant deviation from the SM predictions up to energy scales of several hundred GeV. However, New Physics 
^ ' (NP) models: various left-right symmetric models, extended gauge theories including grand unification theories, 
CO models of composite gauge bosons [l| , some extra dimension scenarios [2] , extra neutral gauge bosons [3] and fermion 

■ compositeness models Q predict various deviations out of the SM predictions and, therefore, their testing in the new 
. energy scale (thousands GeV) is one of the main tasks of modern physics. The forthcoming experiments at the collider 

LHC provides such possibility and probably will shed the light on this important problem in the immediate future. 
J — , The experimental investigation of the continuum for the Drell-Yan production of a dilcpton pair, i.e. data on the 

■ cross section and the forward-backward asymmetry of the reaction 
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pp j,Z ^ fi+^j. X (1) 



O i' at large invariant mass of a dimuon pair (see [3 and references therein) is considered to be one of the powerful tool 
^ in the experiments at the LHC from the NP exploration standpoint. 

The studies of the NP effects are impossible without the exact knowledge of the SM predictions including higher 
order QCD and ElectroWeak radiative Corrections (EWC). The EWC to the reaction ^ are studied well (see papers 
on pure QED corrections jB], and the QED and electroweak corrections in the Z-peak region and above in Q and 
numerous papers cited there). The EWC result contains so-called Double Sudakov Logarithms (DSL) Q, i.e. the 
■ expressions which are growing with the scale of energy, and thus giving one of the main effect in the region of large 
invariant dimuon mass. By now extensive studies have been done in this area. For instance, the weak Sudakov 
expansion for general four-fermion processes has been studied in detail (see, for example, 0] and the recent paper 
along with the extensive list of references therein) . Obviously, the coUinear logarithms of QED radiative corrections 
can compete with the DSL in the investigated region. This important issue has been studied at the one-loop level in 
, where both the QED and weak corrections have been calculated for M < 2 TeV, but it has yet remained unsolved 
in the re gion of M > 2 TeV (with the exception of short numerical estimation of EWC to pp e^e~X by program 
ZGRAD |T]| in see Fig. 7 there). Other important contributions in the investigated reaction at high invariant 
masses are the higher-order corrections (two- loop electroweak logarithms, at least), which also have been studied in 
the works |10|. U M (see also the numerous papers cited therein) . Weak boson emission contribution has been recently 
calculated in [12| and the contribution of higher-order corrections due to multiple photon emission has been computed 
in these contributions are beyond the presented calculations. 

Thus, for the future experiments at LHC aimed at the searches of NP in the reaction ([T]) it is urgent to know exactly 
the SM predictions, including the radiative background, i.e. the processes, which are experimentally indistinguishable 
from ID). The important task is the insertion of this background into the LHC Monte Carlo generators and they 
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FIG. 1: Feynman graphs for the Born (a), one-loop virtual diagrams (b-h) and bremsstrahlung diagrams (i-1). Unsigned helix 
lines mean j or Z. 

should be both accurate and fast. For them to be fast it is necessary to have a set of compact formulas for the EWC. 
They are obtained in our previous paper T5j using the Asymptotic Approach (AA) and improved in subsequent paper 
[TgI . To speed up the bremsstrahlung calculation and increase its accuracy it is very important to have all formulas for 
the cross sections and kinematical restrictions in explicit form appropriate for adaptive multidimensional integration. 
We obtained such form here and reached the high speed and good accuracy of calculation using the 9- and (5-functions 
apparatus. 

II. NOTATIONS AND THE CROSS SECTION WITH THE BORN KINEMATICS 

Our notations are the following (see Fig.l,a): Pi(p2) is the 4-momentum of the first (second) quark or antiquark 
with the flavor q and mass m^; fci(fc2) is the 4-momentum of the final muon fi~^{fj,~) with the mass m; q — ki + k2 
is the 4-momentum of the i-boson with the mass (i = "f, Z); Pa{b) is the 4-momentum of initial nucleon A(B). 
We use the standard set of Mandelstam invariants for the partonic elastic scattering s, t, u: 

s = {pi+ P2)^ t = {pi- kif, u = {ki- P2)^ (2) 

and the invariant S = {Pa + Pb)^ for hadron scattering. The invariant mass of dimuon is M = 

For a start let us present the convolution formula for the total hadronic (H) cross section, where we used such 
abbreviations and indices: "Born" (index 0), V-contribution: (indices: BSE for boson self energies, HV for "heavy" 
vertices, "77" for the IR-finite part of 77-boxes, "7Z" for the IR-finite part of 7Z-boxes, " ZZ" for the ZZ-boxes, 
"WW^^ for the M^W^-boxes, fin" for the sum of "Light" Vertices (LV), the IR-part of 77 and 7Z-boxes and emission 
of "soft" photon, with the energy less than uj). The "/m"-part is IR-finite and described by Born kinematics. Also 
we used common index for V-contribution y = 0, BSE, HV, b, fin and special index for boxes b = 77, jZ, ZZ, WW. 
Thus, the hadronic cross section looks like 

1 1 

a^^^Jdx.J dx2 Jdt [//(^i: (^2, + ft{^i,Q')f,''{^2,Q'H\t)]0{t + s)eM0D, (3) 

here the f^{x, Q^) is the probability at energy scale of constituent q with the fraction x of the hadron's momentum 
in hadron H finding, s = X1X2S, the ^-function under integral sign determined by the kinematics of parton reaction, 
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the factor 

eM = 0{s - Ml)e{Mi - s) (4) 

provides the integration in interval of invariant mass Mi < M < M2 and the factor 

^ e{c ~ cose)0{c + cose)0{c - cosa)e{c + cosa)e{pT{fi+) ~ p^''^)0{PT{f^^) ~ pT'') (5) 

cuts the region of integration according detector geometry, the parameters C* and p™™ will be discussed below. The 
expressions for the angles {9 (a) is the scattering angle of the muon with the 4- momenta ki (/c2) in the center mass 
system of hadrons) and energies (also in the centre of hadron mass system) can be obtained as special situation 
of radiation case ("radiative" invariants v, z, ui, zi=0) from the formulas l|17p presented below. For transverse 
components the expressions take place: _PT(Ai^) = ^losinS, Pt(m") 

— K20 sm a. 

Let us enumerate all quark cross sections in ([3]) using agreement cr(t) = da/dt. 
The Born cross section looks like 

af{t)=^-^ E DW^*{b^t' + l>^:^u% (6) 
where the non-radiative boson propagators look like 



mj + irrijTj 



Tj is the j-boson width, 



7 n,fc \ n.fc \ n.k , \ n.k \ n,k /c>\ 

"± ^ \+ ^1+ (8) 
and the combinations of coupling constants for /-fermion with i- (or j-) boson have the form 

A/^^' = v}v} + a)a), A/^I^' = v)a) + a)v), (9) 

where 



V 



Q f is the electric charge of fermion /, /| is the third component of the weak isospin of fermion /, and sw (cw) is the 

sine(cosine) of the weak mixing angle: sw — \J\ — c^, cw = mw l^z- 
The BSE-part is 

+nfD^ E E i\v^'x' + K-'^^V)Bx]- (11) 

1=7, Z X= + ,- 

Here Ilj''^'''''^ are connected with the renormalized photon-, Z- and 7Z-self energies [13, 13 as 
The HV-part looks like 

Ana 



W^^Re E D^D^* E (A^:'^'x^ +A-Ar-)5„ (12) 

i,j=1,Z X=+,- 



j11 (i\ - 



where the form factors Aj ^"^ are explained in [l 
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The boxes can be presented as 

'^fW = — E D''*{5'^'it,u,b+,b_)-6'^\u,t,b.,b+)), (13) 

where functions 5^'^{t, u, 6+, 6_) and all prescriptions for them can be found in [TBI, [l^. 

The "/^'^''-part (the result of infrared singularity cancellation of 77,7^, LV and "soft" bremsstrahlung) is 

S%, . log i| + ,og -i, - 2 4) + log ^ - 2 4) 

-OgQ,(log J- log - + ^ + log2 - + 4Lz2 — ), 
tu u 6 u u 

Jo = 2{Ql (log - 1) - 2Q,Qz log ^ + (log ^ - l)) , (W) 

where IA2 denotes the Spence dilogarithm. Let us note that correction (5^^^ pgj^ is well known and presented, for 
example, in paper [l^ , and the correction jgp^ can be found in the following way 

To find the cross section for qq-case, it is necessary to change i <-> m in the Born part and QqQi —QqQi in 
"/m"-part. The "hat" in formula ([3]) means only s — > s. 



III. "HARD" PHOTONS 



Let us present the hadronic cross section induced by bremsstrahlung (Fig.l,i-Fig.l,l). Introducing the total phase 
space of reaction as 



1 1 



I^[A] = J dxi J JJJJ dtdvdzdui—^=9{Ru,)0M0D A (15) 







where z = 2kip, v — 2k2P, z\ — 2pip, ui — 2p2P and p is the 4-momentum of a real bremsstrahlung photon. 
The factors d^j and look in "radiative" case slightly different in comparison with "non-radiative" ones: 



61^^ ^e{s^ z~v~ M^)e{M^ -s + z + v), (16) 

and for 9^ we use "non-radiative" expression we should change only the angles and energies: 

/, -, t ^ u + z — Ui 

cosO = 1 H r=, cos a = 1 H 

kiQXiV S k2oXi\'S 

where u ~ v ^ s ~ t and zi = z — ui + v. 

The physical region Q is determined by 6{Ru-^), where R^-^ {Rui is the Gram determinant multiplied by constant 
factor) is described by: 



5 



= f [m^(3s — v) + {s — f )(?Ti^ — s — t + v)] + z[m^(s — v) — rUq^s + v) + st + v{s + t ~ v)], 

Cu^ = z^[{Tn'^ + m'^g-2Tn'^{ml+t-v)~2ml{t + v) + {t-vf] + 

(s - 2t) - m2(2m^ + s + 2t - 2v) + (t - i;)(s + t - v)] + 

+ v^[m^ - 2m^(m^ + s + t - v) + {ml - s - t + vf]. (18) 
Then the total bremsstrahlung cross section have form: 

v3 



X— q—u,d,... i,j—'y,Z 



+ (^2, Q') - X/,^(:Ei, 02)/S(^^ q2)jq^q^^^^^,,- + 



(19) 



Indices I, qk and int mean the origin of emmited photon: lepton, quark and lepton-quark interference, i.e. Final 
State Radiation (FSR), Initial State Radiation (ISR) and their INTerference (INT), correspondingly. The "radiative" 
boson propagators look like 

= —2 — — (20) 

s — z — V — nij + iniji j 

and the expressions R can be found in Appendix A. 

Dissecting the region of integration with the help of function 

2^s 

we divide the cross section in two parts: first one is corresponding to "soft" photons with the energy less then 
uj (it goes to IR singularity cancellation in formula ([14]) of Sect. I) and the second one is corresponding to "hard" 
photons with the energy more then w. We realize the numerical integration of (fTS]) (and, certainly, of ([3])) by Monte 
Carlo routine based on the VEGAS algorithm [20| . 



IV. DISCUSSION OF NUMERICAL RESULTS 



First, we want to demonstrate the independence the results on unphysical parameters: "soft"-"hard" photon 
separator uo (Figl^]) and the quark mass (Fig[3|). In these Figs, we can see the relative corrections 

S^^a^/a^ (21) 

to cross sections integrated over interval of invariant dimuon mass 1 TeV < M < 14 TeV and assuming C* = 1 and 
p-mm _ Q Figl2] shows the w-independence for FSR (left picture), ISR (middle picture) and INT-part (right picture) 
separately in wide range of lo: 10^^ GeV < lo < IQ GeV. We can see also the obvious property for sums the SOFT 
and the HARD parts: |FSR| > |ISR| > |INT|; aU of them are negative. 

For the decision of the quark mass singularity problem we used the MS scheme [2l| and the procedure of linearization 
which is well-grounded in [l^]. After all manipulations the part of cross section which must be subtracted to be 
free from the quark mass dependence has the form (here we used abbreviations Q.S.=QUARK SING, and q{x) = 

1110 

'^Q.S.^\ j dxi j dx2 j dz j dt ^ [[q{xi)/^q{x2, z)9{z - X2) + /^q{xi, z)q{x2)9{z - xi))al'^ [t) + 
-s q=i^,d,--- 



+ {q^q))]e{t + s)eM6D, 



(22) 
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FIG. 3: Independence of the ISR-part of bremsstrahlung cross section on quark mass ruq, the "TOTAL" means "SOFT+HARD- 
QUARK SING.", we suppose here lo = W GeV. 



where 



Ag(x,z) 



2tt 



-q{-,Ml)-q{x,Ml) 
z z 



1 + 

1 - z 



log- 



21og(l - z) - 1 



(23) 



and Msc is the factorization scale [2l| . Figl3] shows the mq-independence for ISR part of cross section (the asterisks 
on plot are the points where the calculation was made, they are connected by straight lines). Wc can see that in the 
range of rather big values (10 - 100) of ratio mg/m„ the difference (SOFT+HARD)- (QUARK SING.) is constant 
(i.e. independent on rUq). In the region of small ruq this property is slightly broken. The reason is simple: at small 
parameter of mass the calculation of mass singularity cross section demands of more time (it is better to say - more 
iterations in integration). In FigOall of points are obtained with the same number of iterations, so in the region of 
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small iriq the result for HARD ISR part has not so good accuracy, in actual fact this part is slightly more. Surely 
increasing the accuracy (and simultaneously the running time of code) we provide the exact cancellation of to^, this 
obvious graph of less importance than Figl3]and we do not present it here. 

In the following using FORTRAN program READY [23| (READY is " Radiative corrEctions to lArge invariant mass 
Drell-Yan process" ) the scale of electroweak radiative corrections and their effect on the obscrvables of the Drell-Yan 
processes for future CMS experiments will be discussed. In FigH] and FigE] it is shown the differential Born cross 
section and the relative corrections to it 

- dM ' dM ^^^> 

as a functions of M. The pure weak (total electroweak) corrections are presented in left (right) picture of FiglS] The 
translation from total to the differential cross sections realized with the help of trick presented in Appendix B. 
We used the following set of prescriptions: 

• investigated reaction is (1) with the energy of LHC = 14 TeV, 

• the set of SM input electroweak parameters: a = 1/137.03599911, mz = 91.1876 GeV, mw = 80.37399 GeV, 
Tz = 2.4924 GeV, Tw = 2.0836 GeV, mn = 115 GeV, 

• muon mass = 0.105658369 GeV, masses of fermions for loop contributions to the BSE: rrie = 0.51099892 keV, 
rrir = 1.77699 GeV, m„ = 0.06983 GeV, = 1.2 GeV, rrit = 174 GeV, ma = 0.06984 GeV, m, = 0.15 GeV, 
rrib = 4.6 GeV (the light quark masses provide Aajfj^(m|)=0.0276), 

• 5 active flavors of quarks in proton, their masses as regulators of the collinear singularity ruq = 10 x m^, 

• non-diagonal elements of CKM matrix = 0, diagonal ones — 1, 

• "soft" -"hard" photon separator u! = 10 GeV, 

• the MRST2004QED set of unpolarized parton distribution functions [l^ (with the choice Q = Mg^ — mz), 

• we impose the experimental restriction conditions on the detected lepton angle — C* ^ C ^ C* the rapidity 
\y{l)\ — 2/(0*1 sse ([5]); for CMS detector the cut values of C* and y{l)* are determined as 

a* 

y{l)* = - log tan — = 2.5, C* = cos 6* w 0.986614, (25) 

also we used the second standard CMS restriction pt{1) > 20 GeV, 

• here we used so-called "bare" setup for muons identification requirements (no smearing, no recombination of 
muon and photon). 

Let us discuss briefly the effects of EWC induced by different contributions (in region M — 1 TeV). The BSE- 
contribution is positive and ^ 0.12, it is usual effect of BSE. The HV-part gives positive contribution ~ 0.07 in 
spite of the negative sigh of DSL (— log'^(m|j^y^/s) = ~^^{w) s) diagrams Fig.l,b and Fig.l,c with the Z and 
W as additional virtual particle. Analysis shows that the Single Sudakov Logs {SSL— I z(w),s) of diagrams Fig.l,d 
and Fig.l,e play the very important role in the region of TeV's M. To determine that we can compare, for example, 
the coefficients at the functions A2{k'^,Mw) (it contains DSL and SSL) and A3{k'^,Mw) (it contains only SSL) in 
formulas (6.8)-(6.12) from 17], the second one is sometimes much more than first one (up to 9 times), whereas 
|DSL/SSL| = |SSL| W4.79 at Af = 1 TeV. The combined effect of all HV becomes positive. Then, the WW boxes are 
uniquely negative and they are the dominant contributions. Let us explain it by the example of WW-diagram and 
7-exchange Born diagram interference: extracting this part of cross section (denote it o'^wx-y) ^^'^ retaining only 
u-type of q uark contributions and leading power of Sudakov logs (there is no SSL in WW-boxes). Then (see formula 
(40) of 

<w% - ■ S'^'^'-'it, u,b+,b^)+uu- ^^^■^(u, t,b+,b^), (26) 

here u(u) = /^(g-)(a;i(2)). In that way we can see the fact: the terms uu and uu contain the same b and different 

invariants t and u as arguments of S^^^'^ . Further, using b^^''^ = -2{v^^ fQq and 6^^'^ = (sec dH]), dH), and 
■y'J = v^v^ + a'a-') we can make sure that 



wwx-y 2(w ) Q^uu ■ t l^y t + UU ■ u l^^ J < 0. (27) 
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FIG. 4: The differential Born cross section as a function of M. 
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FIG. 5: The relative corrections (pure weak corrections - left picture and including the QED corrections - right picture) 
as a functions of M. 



Corresponding contribution of d-type of quarks also less than zero and looks like 



(28) 



The same situation takes place also for WW x Z-ease. At last, ZZ-, ISR-, INT- parts are small enough to give 

determinant effect (ISR gives 0.019, INT gives 0.008, ZZ is ~ 0.0003) and FSR-part is negative and 0.071, 

so the total effect of EWC is found to be negative ~ —0.056. 
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V. CONCLUSIONS 



The complete electroweak radiative 0{a) corrections to the DreU-Yan process at large invariant dimuon mass have 
been studied. For the shortening of code running time (keeping an enough accuracy) we simplify the calculation as 
much as possible (using A A and generalized functions). Using FORTRAN code READY the numerical analysis is 
performed in the high energy region corresponding to the future experiments at the CERN Large Hadron Collider. 
To simulate the detector acceptance we used the standard CMS detector cuts. The radiative corrections arc found to 
become large at high dimuon mass M, the complete corrections at "bare" setup change the dimuon mass distribution 
up to ~ -5.6% (-23.2%; -35.3%) at the LHC energy and M = 1 (3; 5) TeV. 

Some issues have became beyond the scope of the presented paper (the detailed numerical analysis of process 
pp —^ e^e~ X and other interesting observables: total inclusive cross section, forward-backward asymmetries; al last, 
"calo" results - taking into consideration also smearing and recombination). All that will be the subjects of future 
investigation but, first of all, due the importance and complexity of investigated problem, we should cross-check with 
the results of other groups (programs SANC [22], ZGRAD 0,lill, ...) to make sure that our result is correct. Author 
will be grateful to all interested groups for giving a chance to compare the results in that stage. 
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APPENDIX A: EXPRESSIONS FOR THE R 



The expressions for the R have such form: for lepton emission (see Fig.l,i and Fig.lJ) 

2 

= 2s- 2—^{2ul + 2uis + 4uit - 4uiv + + 2st - 2sv + 2f - Atv + 2v'^) 

-2—(-ul - uis - 2uit -s^ - 2st - 2f) - -Uui +4s + 6t- 3v) 
zv z 

-2^(5^ + 2st + 2t'^) - -s(-z + 2ui + 2s + 2t), 

V 



TTt^ S 

Rl"^ = 2^s{2ui + s + 2t - 2v) + 2—{-ui - s - 2t) + -(4s + 2i-t;) 

9 

TTl S 

-\-2—^s{s + 2t) + - {-z + 2ui + 2s + 2t) - 2s, (Al) 



for quark emission (see Fig.l,k and Fig. 1,1) 



9 

Tfl 

Rqkf = 2{z-s + v)- 2—^{z^ + 2zt + + 2st - 2sv + 2t^ - 2tv + v^) 

H 

-2^—i-s^ - 2st + SV- 2t^ + 2tv - u^) 

- — {z'^-zs + 2zt + 2s^ + 2st -SV + At^ - 2tv + v'^) 

Zl 

9 

-2^(z^ - 2zs - 2zt + 2zv + + 2st - 2sv + 2t^ - 2tv + v'^) 

ui 

— -(z'^ - 3zs - 2zt + 2zv + 45^ + 6st - 5sv + At^ - 6tv + 3v^), 

Ul 

i? . 99 = 2—^(-z'^ - 2zt + + 2st - 2sv - 2tv + u^) + 2 (-s -2t + v) 

Zf ZiUi 

+—{-z'^ + ZS- 2zt + 2s^ + &st - isv - 2tv + w^) 

Zl 

2 

777 

+2—^(z'^ - 2zs - 2zt + 2zv + + 2st - 2sv - 2tv + w^) 
ul 

+—(z'^ - 3zs - 2zt + 2zv + As'^ + 6st - 5sv - 2tv + v'^), (A2) 



for lepton-quark interference 



99 = 2{z-Ui-s-At + 3v) + —{2s^ + Ast - 2sv + At^ - 2tv + v"^) 

ZZi 

+ — (2s^ + Qs^t - 6s\ + Sst^ - lAstv + 7sv'^ + At^ - lOt'^v + 9to^ - 3?;^) 

ZUl 

2 

+ -(uis — uiv + + 2st — 3sv — 3tv + 2v'^) 
z 

+ — {z'^s + zH + 2zst + 2zt'^ + 2s^ + 6s^i + Sst^ + At^) + —{zt- 2s^ - Ast + sv + tv) 

ZlV Zl 
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t 2 

H (z^ - 2zs - 2zt + 2s^ + Ast + U^) + -{-zt - uis - - 2st) 

VU\ V 

+ — (-z^ + izs + 5zt - 3zv - 4s^ - 12st + 8sv - 12t'^ + 13tv - 5w^), 

Ul 

R^^^vi ^ 2{-z + 2s -v) —{2s^ + Ast - 2sv - 2tv + v^) 

ZZl 

3_(2s3 + Qs^t - Qs^v + 4st^ - IQstv + Ssw^ - 2i^u + itv'^ - v^) 

ZUi 

- — (s -v) —{~z^s - z^t - 2zst - 2zf + 2s^ + &s^t + Ast^) 

Z ZiV 

— -{-zt - 2s^ - Ast + sv- Af + tv) —{z'^- 2zs - 2zt + 2s^ + Ast) 

Z\ VUl 

2s 1 

(z - s) (~z^ + 3zs + 5zt - 3zv - As^ - 12st + 8sv - At^ + Itv - 'iv^). (A3) 

V Ul 

APPENDIX B: TRANSLATION FROM TOTAL TO THE DIFFERENTIAL CROSS SECTION 

To translate the total non-radiative ^ (or radiative (fT9|) ) cross section to the differential one we should just 
differentiate it on variable M using obvious rule 



dcr da 



dM dMi 
After that and taking into consideration the formula 

de{x) 



(Bl) 

Afi = M 



da; 



5{x) 



we can significantly simplify the form of cross section, as we are in a position to integrate analytically over one of 
variable (here we choose X2) in such a way 

6 

f{x)5{x - z)dx = f{z)9{z - a)e{b - z). 

Finally, we get very simple recipe for translation from total to the differential cross section: so, for radiative case to 
pass 

da" 

H "-"hard 
hard 

we should in formula (fT9|) 

1. do not integrate over X2 and omit 9{s — z — w — Ml), 

2. substitute X2 (M^ + z + v){Sxi)-^ (or s ^ A'P + z + v), 

3. multiply by factor 2M{Sxi)-'^e{Sxi - AP - z - v). 

For nonradiative case the translation steps are the same but "radiative" invariants should equal to zeros: v — z — 
since in that case p 0. 



